Abstract

Precipitation and river discharge are critical components that define hydrological regimes of watersheds. Recently, this region has experienced drought conditions and water resource managers need detailed information about changing flow conditions to better plan for decreased streamflow. The purpose of this report is to provide an analysis of the long term trends in precipitation and discharge in the Breitenbush River watershed in the Oregon Cascades, USA. In this analysis, daily streamflow data will be used to detect changes in the streamflow regime by calculating the annual daily mean for multiple periods from 1932 to 20202. The change in mean will be expressed as percent, slope (%/year), and standard deviation of discharge. R tools from the US Geological Survey will be used to import data and perform the long-term analysis. These tools include dataRetrivial and EGRET.

Introduction

Declining precipitation as snowfall, is causing runoff to occur earlier and more rapidly (Feng and Hu, 2007). This can lead to larger peakflow due to flashy rainfall events. In addition, low flow conditions can be effected by changes in precipitation phase change because snowpack releases water gradually into summer. This can can impact the water resiliency of communities and ecosystems that rely on late season base flow (Rumsey et al., 2015). The Oregon Cascades are a relatively low elevation mountain range and is vulnerable to changes in surface temperature and snow pack reduction. The Breitenbush River subwatershed (HUC 17090005) is located in the North Santiam Watershed in Marion County, Oregon. The climate is warm-summer mediterranean and receives about 59 inches of rain annually. The subwatershed also receives about 146 inches of snowfall annually. Most of the precipitation falls during the winter months. Streamflow monitoring has been ongoing since 1932 and a USGS stream gaging station (44.75262), -122.129) is the primary data source for this study. The upstream drainage area is 279.7 square miles.

Figure 1. Location of Breitenbush subwatershed with stream gage location and upstream drainage area.

Data and Methods

Packages developed by the U.S. Geological Survey were used to determine long-term trends in streamflow. The first is dataRetrieval package (De Cicco et al., 2018) which aids in the reproducability of the results. This package enables rapid acquisition of hydrological data from the National Water Information System which is database containing historical and instantaneous time-series. The entire daily discharge record was retrieved; the use of daily discharge values for computation of historical streamflow statistics is common practice in hydrology. Due to a large data gap from 1987 to 1998 when the gaging station was not active, as a result two periods were analyzed using the EGRET package (Hirsch & Cicco, 2015). This package allows for exploratory data analysis of trends in river data, including water quality trends.

Two aspects of stream flow trends were investigated. The first is the occurrence and variability of 7-day low flow which is the average discharge during the 7 days with the lowest flow for the period specified. Low flow trends are correlated with anthroprogenic water use and precipitation variability (Hammon and Flemming, 2021). Second, maximum daily trends were analyzed during the winter months (November through February).These streamflow metrics were used to assess streamflow response to changes in the precipitation regime.

Results

Maximum daily flows increased substantially through from the beginning of the streamflow record until the late 1960’s when this trend flat lined and eventually began a gradual negative trend that continues to the present day. A similar, but more defined trend is present when 7-day minimum flow is investigated. In this case, an upward trend of 24 percent is observed until the late 1960’s when this trend reverses. A negative 5.5 percent decline in 7-day minimum flow is observed from about 1968 to 1987 and a negative 15 percent decline is record between 2000 and 2019.

1936 to 1987

Figure 2. Plot of maximum daily streamflow metrics for 1936 to 1987.

Figure 3. Plot of 7-day minimum flow from 1934 to 1987.


   BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
   Season Consisting of Jul Aug Sep
    maximum day 

             Streamflow Trends
   time span          change        slope       change        slope
                        cms         cms/yr        %            %/yr
 1932  to  1950          1.5        0.085           14          0.8
 1932  to  1968          2.6        0.072           24         0.68
 1932  to  1987          1.9        0.034           17         0.32
 1950  to  1968          1.1        0.059          8.8         0.49
 1950  to  1987         0.33        0.009          2.7        0.074
 1968  to  1987        -0.74       -0.039         -5.5        -0.29

Table 1. Trend analysis for maximum daily streamflow, 1932 to 1987.


   BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
   Season Consisting of Jul Aug Sep
    7-day minimum 

             Streamflow Trends
   time span          change        slope       change        slope
                        cms         cms/yr        %            %/yr
 1932  to  1950         0.43        0.024           14         0.76
 1932  to  1968         0.59        0.016           19         0.52
 1932  to  1987         0.37       0.0068           12         0.21
 1950  to  1968         0.16       0.0088          4.4         0.24
 1950  to  1987        -0.06      -0.0016         -1.7       -0.045
 1968  to  1987        -0.22       -0.012         -5.8        -0.31

Table 2. Trend analysis for 7-day minimum streamflow, 1932 to 1987.

1998 to 2019

Figure 4. Plot of maximum daily streamflow, 2000 to 2019.

Figure 5. Plot of 7-day minimum streamflow, 2000 to 2020.


   BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
   Season Consisting of Nov Dec Jan Feb
    maximum day 

             Streamflow Trends
   time span          change        slope       change        slope
                        cms         cms/yr        %            %/yr
 2000  to  2010          1.1         0.11            1          0.1
 2000  to  2019         -1.8       -0.093         -1.6       -0.082
 2010  to  2019         -2.9        -0.32         -2.6        -0.28

Table 3. Maximum daily streamflow, 2000 to 2019.


   BREITENBUSH R ABV FRENCH CR NR DETROIT, OR.
   Season Consisting of Jul Aug Sep
    7-day minimum 

             Streamflow Trends
   time span          change        slope       change        slope
                        cms         cms/yr        %            %/yr
 2000  to  2010        -0.28       -0.028         -7.9        -0.79
 2000  to  2019        -0.54       -0.028          -15         -0.8
 2010  to  2019        -0.26       -0.029         -7.9        -0.87

Table 4. 7-day minimum streamflow, 2000 to 2019.

Conclusion

The Breitenbush River has seen periods of increasing streamflow metrics during the first half of the 20th century, but recent trends indicate a reversal of these trends. Water resource managers downstream of this basin should be aware of the severity of this negative trend, because it may have implications for water security. Communities that rely on summer time low flows for irrigation could faces reduced water resource resiliency. As Oregon enters another year of record drought, it is important to note the historical trends do not indicate that we will return to normal conditions in the near future.

Work Cited

De Cicco, L. A., Lorenz, D., Hirsch, R. M., & Watkins, W. (2018). dataRetrieval: R packages for discovering and retrieving water data available from US federal hydrologic web services. Reston, VA, DOI, 10, P9X4L3GE.

Feng, S., & Hu, Q. (2007). Changes in winter snowfall/precipitation ratio in the contiguous United States. Journal of Geophysical Research: Atmospheres, 112(D15).

Hirsch, R. M., & De Cicco, L. A. (2015). User guide to Exploration and Graphics for RivEr Trends (EGRET) and dataRetrieval: R packages for hydrologic data (No. 4-A10). US Geological Survey.

Rumsey, C. A., Miller, M. P., Susong, D. D., Tillman, F. D., & Anning, D. W. (2015). Regional scale estimates of baseflow and factors influencing baseflow in the Upper Colorado River Basin. Journal of Hydrology: Regional Studies, 4, 91-107.